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1 . INTRODUCTION 

Applied mathematics often involves*" the solution 
of single' equations or systems of equations. In this 
nodule we study a particular technique for sofving 
equations in the context of functions on the real line,. 
The technique includes an iterative'process. i.e. , an 
algorithm, that can be understood and interpreted with 
>ease pictor*ally and can be* implemented readily on 
inexpensive and widely availaole calculators (program- 
ming feature helpful but not necessary). We also 
present some examples of ^physical problems in which 
these, equations arise, « 

. ' 2. THE PROBLEM 

The basic problem here is to solve equations" of 
"the* form ■ \ c 

CD f(x) « x, * 

A solution of Equation (1) is called a fixed point of 
f|^ geometrically, such a solution is a value of x for 
which the graph of the function y * f (x) crosses the 
line y * x (see Figure 1), For example, if f(x) = x 2 , 




\ Figure 1, Geometric interpretation 

of a fixed point* 



then Equation (1) takes the form 



x 2 ■ x, ■ 

which has two solutions (fixed points), x = 0 and 

x = 1, . - 

Equation (1) is more general than it may seem* at 
first, because' an? equation of th*e* f orm y 

(2) g(x) * 0 

may be converted to the form (1) by adding x to both - 
sides: 

ti) g(x) + x = x, 

and Equation (3), which is equivalent to (2), is now 
of the form (1), with f/x) = g(x) + x, 7 The method 
presented here may therefore be applied to find roots 
of equations, ay well as fixed points, * 

V THE ALGORITHM ' - 

1 

The iterative process that constitutes the basis of 
the contraction mapping principle is often called Heard* 6 
algorithm or simply iteration* The process is easy to 
describe: for a function f whose range is contained in 
its domain, start with a point xo in the domain, and then 
successively reapply the function. That. is, let 
* x n*n=l be the sec l ue nce i % 

(4) x x * f(x o ), x 2 = f( Xi ), x 3 = f(x z ), y ,. 

Under some/simple conditions that involve the con- 
cept of a contraction (to be described later), the sequence 
defined in (4) converges to a fixed* point of f(x); that 
is, if we let x = lim x , then 

f<x) = x. 

We shall formulate the definition of a corTtTaxrron— — 
and a statement of contraction mapping principle in Sec- 
tion 4, but first we illustrate the algorithm by some 
simple' examples , 

% 



~ , - 3,1, Example 1 ' >* 

. \ ' " • ». - . * * 

. Let irf begin .with an example that is easy to see 

and for which we know the outcome in advance. Let 

f(x) » jx + 2, and« suppose we wish to solve the equation 

K-' - ■: . + - x . ..." . » . . . . 

(see figure. 2): By elementary algebra, we know that r ~~" 

I, 



* i 




^ X 



Figure 2, Geometrical interpretation 
of Example 1 , 



the solution is x = 3. To apply the PicarS algorithm 
we choose any starting point, say x 0 = -5, and calculate 
successive values of the function according to,>{4),' 
The results are indicated in Table 1; check them on 
your calculator.. 

. -The fixed poiijt x = 3 is reached correct to four 
. placds in eleven steps. (Note that had we rounded to, two 
.places, -we would have reached the solution to two places in 
only five steps;)* 
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TABLE 1 

Results of the Picard Process* Appi ieJ to 
the Function in Example 1, wjth x 0 * -5, 



» 


V 

t 

Formula for f fx) 


Numerical / 
Approximation 
of ffxi . 


-x 0 ~ = -5 


i U5) + 2 ' 










i (.3333) + 2 * 


9 1111 bp • 
£ *±±±± = X2 


J 2 


▼ 12,1111) + 2i 


Z , / {Jo / 1 = x 3 




i f 2. 7037 ) + 2 V 


O Qni -> 




4 C2. 9012V + 2 * x a 


O Q*C 7 1 

f 




j (2,9671) + 2 


2,9890 , 1 




j (2.9890) + 2 


2,9963 




y (2,99^3) + 2 


2,9988 




y (2,9988) + 2 


2,9960 


I 


,y (2,9960) + 2 


2,9999 


x 10 


*y* (2,9999) + |2 


3,0000 



It is ^important to understand the geometry of the 
Picardj procedure* When we let x l « f(x 0 ) and then calcu- 
late f(x x ), the first y-value becomes the new x-value, / 
and then we move tt> a new y-value. This process can be 
represented by a path in which we start at (x 0 , 0), 
move vertically to (x 0 f f(x 9 ))' or (i Q , x 2 ), then move > . 
horizontally to (x Jf x^', and then vertically to 

£(x a )) f and so on, A picture of this procedure 
is' shown in Figure 3,- 




4 



•7 



y ' 



3 _<Vf<x,)) 



(x 0 .f(x 0 ))- (x 0 ,x,) 




2. till 



Figure 3. Centric Interpretation of the P.card' procedure for 
the function In Example 1, with x Q - -5. 



Exercise 1/ 

Use the P?card procedure to solVe? the equation ic - 1 . x . 



3.2 Example* 2 ' , ' 

I .'Suppose we wish to solve the 'equation 

COS X » X. / 

In this case, the equation, is transcendental, and 

algebraic manipulations do not yield a solution. . 
/•A.reasonable approximation to the solution can be 

obtained graphically. Or, if you have tables for 
. trigonometric functions handy/ you may simply scan 

for the poin't where cos x - x (radians) ' The CRC 
, tables shoi, that cos(.73) i ..7452 and cos(.74) i .7385 
■ so we would expect the solution to be close to .74. 

The Picard procedure for the equation cos x = x 
•can be implemented easily on many calculators. With 
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,t$e ^calculator in the radian mode, start with x = 0, 
and repeatedly calculate the cosine (see Table 2). 

* * • 

TABLE 2 , 
Nuiaerlcal Solution of the'fEquation cos'x^ x. 



X 0 = 0 



23 



Formula for f(x) ' 
cos 

cos (/5403) 

cos (.857 3) 

cos (.6543) 

cos (.7935) 

,cos (.7034) < 

cos (.7640) 

cos (.7221) 

cos (.750*4) , * 



cos (.7391) 



= x, 



Numerical 
Approximation 

- of f-(-x4 

f.0000 
.5403 V 
.8576 
. .6543 
'.7935 
.7014 
.7640 
.7221 
.7504 
.7314 



,7393 



The convergence for the equation of Exampie 2 was 
much slower than for Example 1; this time 23 iterations 
were required for four place accuracy (11 for two place 
accuracy) . You may be able to see the reason for the 
slox convergence by studying the geometric picture 
here we "wind^around" the fixed point, as shown in 
Figure 4% 

\ > 




■. Figure k 



Graph showing the Hicard procure* for the 
equation cos x ■ x. * 



3.3 Example 3 0 ' " * 

Consider a. rigid bar that is hinged at lane end, 
with the other end free to swing. The bar has length 
* and, uniform mass distribution, with total mass'm, 
Let f 6 8(t) bethe ang-le the bar makes wit^j the vertical 
Vat time t. At the hinged end^of the bar, there is a 

.torsion spring that exerts -a force> of magnitude ke , where 
-«j^4A_a^nstA^ vertical ^rest^ 

position. .We suppose also that St downward .IfarcejJE acts 
at th0 free end of the bar* The* bar is given an initial 



12 




Inge, and spVlng 



Figure 5. Illustration for 
Example 3* 



velocity 6 f and started into motion at time t = 6, 
o 

(See Figure Z.j 

If we neglect the gravitational force, the 
equation of motion can be obtained by considering 
(angular) forces; this equation is 

(5) jm* 2 8 - k8 + JtF sin 6 « 0, . 

(The factor jmZ 2 is the moment of inertia of the bar 
about the hinged end,) If k < xVF, then the system has 
exactly one equilibrium position between 0 and it. 
By Equation (5), with 9 set equal to 2ero, the equili- 
brium position is characterized as the solution 
(between 0 and tt) of the equation t 

-k8 + £F sin 9*0, 



Thus, the value of 9 that yields equilibrium is the , 
positive solution of the equation 



*F 
X 


sin x * x. 




For example, 


if the values of i, F and 


k are such that 


JtF/k = 1,2, 


then we can find the equilibrium position 


bV the Picard algorithm, starting at x 






TABLE 3 






Numerical Results of the Ptcard Sequence 
for 1 .2 sin" x « x, Starting at x Q « tt/6. * 


x 


Formula for f(x) 


Numerical Approx- 
imation of f(x) 

v o?6ooa = Xj 


X 0 « 7T/6 


1.2 sin (tt/6) 




1,2 sin (.6). 


.6776 = ^ 


x 2 


1,2 sin (,6776) 


,7523 




- 1,2 sin (.7523) 


.8200 


♦ \ 


. mt 1,2 sin (,8200) 


.J77Z 


1 

22 


1.2 sin (1.0267) - 


/ 1/0267 



Thus, in this case, the equilibrium position occurs near 
1.0267 radians, or around 59° ♦ 



Exercise 2 

Find the equilibrium position for the system in Example 3 
When if/k • 1.4. 



4. CONTRACTION MAPPINGS 

In this section we consider conditions under which 
the Picard algorithm can be used to solve" equations* 

' 4*3 Example 4 

Consider the function f(x) «~2x' +1, We know by 
simple algebra that x » i is the solution of the 
equation 2x - 1 » x. However, when we try the algorithm 
with x 0 > 1 or with x 0 < 1, we move away from the solu- 
tion (see Figure 6, and compare Example 1 and Figure 3), 



y 




9 



It turns our that this behavior is related to the fact 
that the slope of the line y = |x + 2 is less than 1 
and the slope of the line y = 2x - 1 is greater than 1, 
as we see in the discussion below. 

Let f be a function that is defined on an interval 
I which is of any ot the following types: 

• {x: a < x <W, where a < b are finite numbers; 

• {x: x < b}, where b is finite; *> 

• {x: x >^ a}, where a is finite; 

• the entire real line, * 

Thus, we consider the interval I as the domain of f 
(even though it may be that the formula for f can be 
applied to a largej set). When the values f(x) are 
also contained in the interval r for x in I, we say • 
that "f maps I into I" and write f: N I, (For 
example, let f(x) = |x 2 , and consider I = {x: 0 < x < 1}. 
If x e I, then f(x) e I also, and so f: I + I, i n 
this example, we note that the natural domain of f is 
the whole Une, but we simply want to consider the 
behavior of f(x) for values of x from 0 to 1.) 

Let I be a real interval as described above, and 
let f: I I. Then f is called a contraction on I 
if there is a constant r such that 0 < r < 1 and for 
all x,t in I we have - 

|f(x) - f(t)| < r | x - t |. v 

4.2 Example 5 . 

Let f(x) = ^x + 2, I the\whole line. Then for 
all x,t we have 

I'M - f(t)| - |[jx ♦ 2) - (ft ♦ 2)| = ||x - t|. . 
Thus, f is a contraction on I, with r = |. 



10 



' Let f (x) - 2x - 1, I the whole line. Then for 

• "any x,t we have 

* i 

Jf(x) - ; f(t)| « |(2x-.l) . (2t-l)j - 2|x-t 
'Thus, f is not a contraction- on I, 

The- contraction mapping principle for the 
real line* Let I be a real interval of " a type 
indicated above, and f: I I, If f a con- 
traction on t, then the equation 

- — r - -~\ f(x) = x 

has a unique solution in I, and this solution 
may be obtained by choosing any point x Q in I, 
forming the sequence 

" f (* e )> x 2 * f (*iJt * 3 ■ 
ana passing to the limit* The solution is 

lim x . 

.. . . n-*°° • 

A proof of the contraction mapping :principle-.is 
.preseated in the Appendix* 

*\ .The contraction mapping principle may~seem to be 
restrictive, but there are two important techniques 
for extending, the applicability. First, we note that 
if f is defined on I, ifVf is an 'expansion 9 i.e. for* 
all x,t we have 

l*M - i(t)[ > k|x - t|\ . , 

wtfere k > 1, ahd if "the range R of f contains I, then 
we must have the following: 

"f is one-to one, so that, f* 1 is defined on R; 

• - f~*: R + I, and I e R, so f ~ l : R R; 

• f" 1 is a contraction on R, with r « 1/k, 



16 



Hence there must be a unique x in R such that 
'f" X (x) = x. 

But since the graphs of f and f" 1 are symmetric in the 
line y x, x must also-be in I and we must have 

f(x) = 
4.4 Example 7. 

Let f(x) « 2x - 1, I the whole line. We saw in 
Example 6 that f is an expansion, with k = 2. The 
inverse function is defined also on I, and is given 
by the formula 

f" 1 (x) - \x + 

Then f 1 is a contraction on I, with constant 
r * 1 = ]c* * The p icard algorithm now yields the 
solution of the equation 

.fx + \ « x, . 

which is-x ■ 1, as shown in Figure 7. (Compare and 



y- 





, * / 

/*- x * y 










x o < 1 


1/ i i 
2—1 1— 


i I, 

^ / 
/ 

/ 

/ 


' \. Xfl>1 



Figure 7. The PIcard algorithm for Example 7. 



contrast with Figure 6.)* Since the equation 
nay be rewritten 

- 2x - i » x;- 

we have the solution, of the original probleirtujosed. 
in Srtalhple 4., , 

The second ,way in which we can extend the Applica- 
bility of the Picard algorithm \is by considering the 
iterates of f rather than £ itself. For f : I + -I we 
let 

• - f/W - f(x), f z (x) - fff/x)), fji) = f(f s (x)).,... 

It turns out that if there exists a positive integer n 
— such-tliAt - £ n - is a contraction on I, then the fixed point 
°C f n is also the fixed point of f in I, and the fixed 
point can be found, as before, by choosing any starting^ 
point x 0 and iterating t;he function f, even though f 
itself is not. a contraction. (It would be a challenging 
exercise for an advanced student who can handle the proof 
— of-^he contraction mapping principle in the Appendix to 
construct a proof for this extension of the principle.) 
; This technique may also be combined with the use of in- 
verses, \as- the following example will illustrate, 
v , 

4.5 Example 8 

-. * * " ~ ■ 

Suppose we wish to solve the equation 

A rough sketch shows that we should expect to" find two 
solutions, one between 0 and 1, the other greater than 
two. x In Figure 8 we can s$e that if we take x 0 to be, any 
point to the left of .the larger solution, the Picard " 

sequence-vri-l-converge -to-the-smallei-solution.- -The 

approximate values obtained by starting at x 0 8 0 are: 
0, .25, .3210, .3446, .3529, '.3558, .3568, .3572, .3573, 
.3574. 

. * i3 



The reason that the-above^ sequence "converges can be 
can be seen by choosing our interval I carefully and 
observing that f is a~ contraction on I. For example, 



y 




Figure 8. The graph for Example 8. 

if we let 

I = {x: x < 1} , 
then f: I + I, and for x in I we have 

f'W =Je X . , - 
Thus, by the law of the mean, for x, t < 1, we have 

" If 00 - f(t)| = |(i e e) (x - t)|, \ 

where £ is between x and t. Thus, 

|f(x) - f(t)| = }e ? |x - t| '< f |x - t|. 

THus f is a contraction on I with contraction constant 
r * e/4. - 

If we^ try to re ach the larger solution of je x = x 
by taking x 0 between tire~solutibns ,~welTustrf aff, T as 
shown also in Figure 8. * This behavior can also be 
understood in terms of contraction mappings. This 

10 



v time we* take as an example* I \« {x: x < 2}. If we 
iterate-f, letting $;(x)" » f(x), f 2 ( x ) - f(f;(x)), 
, ■ f if£M) >♦ • ♦ t as discussed"" above, then it turns 

out that f k is a contraction on I, and in fact for 
x, t < 2,, it can be shown that 

1**00 - f*(t)| < .9|x - t|. 

* • 
• Thus, there is a unique fixed point for |e x in "I and, 
since we have already found one (approximately) , the Picard 
algorithm cannot yield the £ixed % point to the right of^x = 2, 

In^ addition, if we try to find the second solution of 
the equation^ e x = x by choosing x 0 to the right of this 
• | root (try it for yourself, with x 0 = 3, for exampl^) , the 
sequence of iterates diverges, increasing without bound 
exponentially. This behavior,' too, is indicated in Figure 8. 

Thus, if we apply the Picard algorithm to any value 
x 0 (other than the , second oot itself exactly!)*, the se- 
quence of iterates will not converge to the second, or 
larger, root of-£e x = x. However, we can find this second 



•a 

y 




X * o 

H Figure 9. The graph of y - Inkx. 



root by applying the Picard algorithm to the inverse func- 
tion, in 4x , as described on pp. 11-^12. For example, if * 
we take I * {x: x * 1.4} (1.4 is just a little larger, than 
*n4, where the slope of |e x is 1), then f is an expansion- 
on I, and the set of values assumed by |e x for x in I con- 
tains I.j Thus, the inverse should be a contraction, and 
is, as can be seen in Figure 9, 

Now if we take I « (x:x > 1.5}, for example, then 
*n4x is a contraction^ on I, and we may solve the equation 
£n 4x = x 

by the Picard algorithm, starting at any point to the 
right of x « 1.5. If we start at x « 2, we obtain the 
sequence 2, 2.0794, 2.1184, 2.1370, 2^1457, 2.1497, 
2.1516, 2.1525, 2.1529, 2.1531, 2.1532, 2.1533, 2.1533, 
.... Thus, 



*n[4(2.1533)] ='2.1533, 



so that 



e 2.1533 j 4(2.1533), 

1 A 2.1533 • , im 

^e « 2.1533. i 

Thus, the larger fixed point of ^e x is approximately 
2.1533/ 

Note .that in choosing' a value for the left endpoint 
of the interval I indicated above, it was not really 
necessary to ensure that ln$x would be a contraction on I. 
We chose 1.5, which is to the right of the point x ■ 1 
where the curve ln4x has slope one, and since tl>e slope is 
decreasing the mean value theorem shows that *n4x is a 
contraction on I for the choice we jnade.- However, we 
could have chosen instead "any value to the right of 0'.3575, 
which is slightly larger than the first root of the equa- 
tion- *n4x = x. (Recall that our approximation for this 
root, which we obtained using the original equation 
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1 X - * 

|e * x, was 0,3574 to. four decimal- places.) For 

example, if we had taken l"~ {x: x > .4}, then the 
iterate f 5 » would be a contraction on I, so the iteration 
procedure would still yield the desired fixed point; 
this example illustrates the combining of the two exten- 
sions of the basic contraction mapping principle, as 
promised in the final sentence immediately preceding 
Example 8, (Try iterating the function *n4x starting with 
^x 0 » ,4, and follow the^iterationboth on the calculator^ 
'and. geometrically , on a graph of the function--you will- see 
the points, coming closer together starting with the fifth' 
application of the function.) 



Exercise 3 

Solve the equation ycoshx * x. (Hint: This problem Is 
rtnrttarr-to-^^ n^y need t0 know that 

for f(x) « jcoshx the appropriate Invers"r?unHlon for this 
problem is given by the formula f~ V (x) « in (2x + Ax a " l) in order 
to flpd the larger solufclqn; On some calculators the Inverse 
Jwerbblic cosine Is available directly.) 



5. 'FURTHFR EXAMPLES AND APPLICATIONS " 

S,l Example 9 f \ \ 

When a beam of light passes through a narrow slit, 
it spreads out in the shadow region. This effect is 
called diffraction, and a diagram representing a simple 
model, called Fraunhofer diffraction is shown in Figure 
10. 
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By an application of basic principles of optics, 
it can be shown that the light intensities on' the screen 
can be ^approximately) expressed in the form 



I - a 2 sin p 



where 6 is a suitably chosen spatial variable. (See 
D. Halliday and R. Resnick, Physics .) The quantity 



A * A sin g 

0 



is called the amplitude of the vibration. A problem of 
-interest in optics is to determine the~4ocation\of the 
maximum intensities, thus we wish to optimize the function 



A sinfi 



9 0 - J 18 

. . • V - 



We "have 



Begs B - sin 6 ' 
B 2 



^ and we find initial points by setting A 1 = 0: 



B 

The critical values ban be obtained by solving the. equation 
^ ,tan B* B, or tanx x. X 



The geometric picture is shown in Figure \l 



y » x 



y ■ tan x 




Figure 11. The graphs of y » tanx and y » x. 



We see the first maximum intensity at x ^ 0. 



For larger Values of'x the fixed points will be near 
pdd .multiples of *?2, such as 7*/2, 9w/2, ♦ We 
find- the location of the first nonzero intensity here, 
k using the Picard algorithm applied to the function 
y * n + arctan x, starting at x 0 = 0 (see Figure 12) , 
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Fixed point 

4 Figure 12. The Picard algorithm for'flndlng the smallest 

positive fixed point of tanx f using the Inverse* 



The resulting sequence of approximation is 0, 3,1*16, 
4/4042, 4,4891, 4.4932, 4.4934 MAM. The location 

of the first maximum intensity to the right 0 corresponds 
to B * 4,4934, 



Exercise k 

Find the fixed point of tanx that lies between 3it/2 and 
5*/2 (i.e;, the location of\the second maximum intensity to the 

right of 0) In the Fraunhofer diffraction pattern. 

• * 

Exercise. 5* 

In cavity resonators used In travel I ng-wave* tubes the energy 
storage ratio wJ/W 2 has an expression of the form 

W * „ 2s In 2 (972) > ft . % 
where 0 t is constant. Use the Picard algorithm to find the 



See H,M. Nodelman and F."W, Smith, Mathematics for Electronics . 
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s^n x 

«axi«u» v«?ue of—— (and then muttfply the critical value of 
x ky H 2) to f ind the value of 8 that maximizes the energy, storage 
ratio* \" 



S.2 Example 1ft /' 

~/fin submarine- location problems i it is often 
important to. find tKe submarine's closest point- of 



- Appr oa ch (CPA) to a ~" S onobuoy in the water I Suppos.e 
that a sonobuoy is located" at (2, -|>) on a rectangular: 
system and that a submarine travels on a parabolic - * 4 - 
path,- along t the curve . y k 2 . " (See Figure" 13. ) For 
any -point (x,x 2 ) on the parabola, -the distance. jto the 
sonobuoy is- 

D « Ux - 2> 2 ♦ (x 2 + £) 2 ]* - 



(x* + 2x 2 - 4x + 4r) 



1M 



We wish to find the critical x that minimizes D; this 



y « x" 
Submarine track 
In two dimensions. 




Sonobuoy 



Figure J3* Diagram for. Example 11, 



value will also ^minimize D 2 , and we proceed 

Sc (D ^' = 4x3 + " 4x " 4 > ^ 
- so we wish to solve the equation * 
x 3 + x - 1 * 0. 
We can arrange this as a fixed point problem: 
x 3 + x - 1 ' \ 

x(x 2 t 1) « 1 

" * 1 

= x. — - 

1 + x2 "\ 

A rough ^sketch of " the graph of f(x) * j * x2 (see 
Figure 14)-~shows that f is a concentration mapping 
on the whole line, so we implement thevPicard algorithm, 
starting at x 0 = 0, to obtain the sequence 0, 1, .5", ■ 

I 



y 




Figure }k. Graph of y * — - — =• . 

1 + x 2 
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• 8, .6098, .7290,", .6530, .7011, .6705, .6899, 
.6775, .6854, .6804, ,6823, .6823,./..'.. . * 

. Thus the CPA of the submarine to the sonobuoy is approxi 
. mately the point (.6823, (.6823) 2 ) or (-.6823, .4656). - 

--^--toe-conclude this-section-with an~example~that is 
presented in Wylpg's. Advanced Engineering Mathematics . 
5.3 Example 11 

< A sleir(fei-TOinjf~ length ~L-" has iis> curv e d - su r fa ce 

perfectly insulated against heat flow. The rod is 
v. located along- the x-axis with its left end at x = 0, 
... the right at x « L. The left end is maintained at 
. temperature. u(0",t) « 0 for all time t > 0, the right 
end radiates heat freely into air of constant tempera- 
^ ture 0. If the initial temperature distribution in 
the rod is given by 

u(x,0) m g( x ), 

the problem is to find the temperature u(x,t) at any 
point 0 < x < L and t > 0. By an application of Stefan's 
law and a "separation .of variables" argument, Wylie * 
shows that* the* solution can be expressed in the form 

/ (6) u(x,t) - ? B n e - A nt/* 2 sinA x .. 

. n*l . n „ - n 

What is important for us is the fact that the values 
, A n are" given by 



where the_z K ^z 2 , are the positive roots of an* 
equation of the form . A 

* tan z * -az 

~ "where a Ts~T~po s ftive "con tantTsee Figure^ 15) . For 
example^ if a « .25, then we must solve the* equation 
of the form 




Figure 15. Sketch showing z }9 z v z^ for Example 12. 



.25 



tan x 



in order to find the eigenvalues A , A 2 , .... We 
obtain the Picard sequence from the algorithm . 

y 3 it + arctan (-.25x): 

0, 3.1416, 2U758, 2.5873, 2. 5674, 2.5710, 2.5703, 

2.5704, 2.5704, Thus z l i 2.5704, and the first 

eigenvalue can be obtained from the relation 



2.5704 



When we have found several values of A- , we 

n 

may substitute them into Equation (6) and thereby 
find series approximations. for._the- temperature 
distribution u(x,t). ^ 



Exercise 6 

Find the second root z 2 In Example H with a * .25. 



6. COMPARISON WITH NEWTON'S METHOD 

* * 

There are*many numerical methods of salving real . 
equations of "the form 

g(x) - 0. 

(This is equivalent to the'problem of solving . 
f(x) - x, 

since we:may writethe latter as - 

, f cx) - * o : 

and dimply let g(x) - f(x) - x.) One popular such 
method is called Neuton'e Method, which is carried out 
by choosing a starting point*x 0 carefully and then 
applying the algorithm * v 

f ( x n> ■ 

x n+i " x n " rrs^i * 

Newton's method enjoys many advantages; for example 
the resulting numerical sequence often converges faster 
than that obtained from the contraction mapping princi- 
ple. There- is rather general agreement that Newton's 
method is the superior of the two. (A comparison of 
these methods can be fouijd in M.L. Dertouzos, et al,,* 
Systems. Networks, and Computation: " Basic Concepts, ) 

If Newton's method is better-, you may well ask, 
"Why bother with Heard at all?" There are two primary 
reasons., First, there is really not a dichotomy between 
these methods, i.e. you need not choose between them* 
Learn both!' Because it is simpler and -often, easier 
to implement. on a (programmable or nonprogrammable) 

-jgjg^glgg >. the^c. ontga^ionjapping ,principle.serves 

as an excellent introduction to numerical methods. 
The simplicity and ease of implementation, coupled 
with its range of applicability make the Picard 

'* : " v ; ' - 

' - r 30 \ , • \ . .. 



algorithm also a viable alternative to Newton's method, 
in many problems. 

The second reason for learning this method is 
that the contraction mapping principle. can be formu- 
lated in more general terms to solve a. much wider class 
of problems than real; equations. For example, it can be 
used to solve systems of linear equatioiis, nonlinear 
equations in higher dimensions, differential equations, 
and integral equations, (See, for example A, W, Naylor 

Engineering on differential- equations ; they also have ; 
a nice discussion of the relation of the principle to 
closed loop feedback systems,) Concrete applications on 
the real line serve as a solid foundatio^ for exploring 
interesting applications in more advanced settings. 



\ • 
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7. MODEL EXAM 

Solve the following equations uring the contraction mapping 
•prlwcrfphrtP ftmPir a 1 go r i thm) : " 



a, * e~ x - x; 



c. sin.2x 



Use the contraction mapping principle, ♦.ogether with a little 



a. 



b. cot x » x (find the solution that lies between 0 
and it). \ 



manipulation, to solve the equations: \ 



\ • 
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8. ANSWERS TO EXERCISES 

-2.000<h ' x 
1.3720. 

.589<U. 2.1268. 
.7253. 
2.3311. 
5.35*0. 



\ 9. SOLUTIONS FOR MODEL EXAM 

a. .5671; 

b. -7531; 

c. .9477* 

a. Use ^~ - x; then x - .8767. 

b. Use arctan — * x; then a * .8603. 



APPENDIX : 



; L t : : r . : : &*P*£ ^ th^ ^aaip^dont^aoH^n napping principle. 
7 - ;^ that-t Is m interval of the 

^^&r» {x!;a ^ * '< t&fli '-ft X-* I, ind that there 
exists a positive constant, such that for all x, t 
in I Ve have 



r £(t)j < rjoc - tj, 



Let"x 0 be any point in X, and let ~^ 

.*i *(x 0 ),* x a f(X t ), x 5 « f(x a ) .... 

^Tfien7~~~~ ----- 

: l*, : * r lj" M(x;)-f(x,)| < r| Xj -x|j, 

I*, -x 2 | « |f(x 2 )^f( Xl )| < r|x t -^|<r»|x 1 .x t |, 

- J* % -x t | * |f(x 3 ).f(x 2 )| <Jr|x $ -xJ<rMx 1 -xJ i 

It can be shown by mathematical induction that for 
every n, 

Thus, for all positive integers n and k we have . 

- 'J lx n+k - xj < |x n+1 - xj + ... + |x n+k - Wl | 

* < (r n + ... ♦ r"^" 1 ) | ?1 - xj . • 

n+k-1 . ' 
< I?, " x 0 h I r 1 
i»n 

where tS B }* M0 is the sequence of partial sums for the 
series ^' 

oo 

- k fl • 
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CerIc 



"Since this sequence is a Cauchy sequence (remember 
0 < r < l; so the, geometric scries converges), it 
follows that ths sequence . 



^\n*n»l 



is also Cauchy,, and hence converges. Let 



lira x„ ■ "x. 



— ^ow^he-condi^on^hat-^be-a^ontractdT)n-al^o ' 

implies that' f is continuous (for a given e > 0 take 
$ » e in applying the definition of continuity). 
Therefore, 



f(x) = f (lim x ) lim f (x n ) = lin x n+1 = x, 

so x is a fixed point of f . Moreover the fixed point 
of f (in I) is unique, because if x 1 is any fixed point 
of f in I, then * 

Ix - x»| = |f(x) - f(x»)| < r|x - x»|. 
from which 

(1 - |r) |x - x» | < 0. 
But 1 - r > 0, and since |x - x^^O, we must have 
. 0 < (1 - r)|x - x'| < 0, ' 
from which 

(1 - r)|x - x'| » 0 
and finally » 

|x - x» | = 0 ' | 

so that x = x», and the fixed point is unique. 
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STUDENT FORM 1 
Request for Help 



Return to: 

EDC/UMAP 

55 Chapel, St. 

Newton, MA 02160 



Studen^: If you .have trouble with a specific part of this unit olease fin 

yo U u Till wm hefo^ 6 * £° y ° Ur inStrUCt ° r f ° r ^istance?1 n ; £ 
you give will help the author to revise the unit. 



Your Name 



Page 



O Upper 
OMiddle 
O Lower 



OR 



Section 



Paragraph^ 



OR 



Description of Difficulty: (Please be specific) 



Unit No. 



Model Exam 
Problem No, 

Text 
Problem No. 



r 



Instructor: Please indicate your resolution pf the difficulty in this box. 
^) Corrected errors in materials. List corrections "here: 



o 



Gave student better explanation, example, or procedure than in unit, 
" Give brief outline of your addition here: 



Assisted student in 'acquiring general learning and problem-solving 
skills (not using examples from this unit.) 



OP 



9 

.ERIC 



Instructor's Signature 



Please use reverse if necessary. 



STUDENT FORM 2 ' 

• ° • Unit Questionnaire 5 ^ Cna Pel St. 

Naae Newton, MA 02160 

*** ■ s , — ' : Unit No. pate. " 

Institution,^ C onrsP NO. 

Check the. choice for each question that comes closest to your personal opinion. 
*• How useful was the amou nt of detail in the unit? f 



_Not enough detail to understand the unit 
.Unit would have been clearer" with wore detail 
.Appropriate amount of detail i 

-?oi t «T.W' J"*;^ 11 ^ to ° ^tailed, but' this was not distracting' 
Joo much detail; I was often distracted 



2 \ - How helpful were the. problem answers? * 

■fiK.S^SSLSEL'SL'S!?* 1 co »" «*•*• i— step. 



_ — S^f f icient information was .given to solve the problems 
Sample solutions were too detailed; I didn't need them 

3 * * ^L f0 !J^ f ^ llng /? P ^ rerequisltefl - >" >" much did vou use other source, (for 
exam ple r instructor, friends, or other h o oks) in oriir to understand thj mUrt - 

A Lot Somewhat _A Little J Not at all 

A# l°lJ£«l n S !: hlS Un i\ in com P arlson t0 the amount of time you generally spend on 
a lesion (lecture and homework assignment) in a typical math or science course? 

Jfuch Somewhat About Somewhat Much 
Longer Longer . the Sme shorter Shorter 

5 * Were any of the following parts of th e unit confusing or distracti ng? (Check 
as many as applyT) • — ' ft 

* Prerequisites 

- Statement of skills and concepts (objectives) 

Paragraph headings 

. Examples / 

Special Assistance Supplement (if present) 

„ Other, please explain " *" 



6 * gere any of the following parts of the unit particularly helpful? (Check as many 

as apply*) ~ * J 

Prerequisites 

Statement of skills and concepts (objectives) 

Examples 

Problems » 

_ Paragraph headings 

Table of Contents 

Special Assistance Supplement (if present) / 

• Other, please explain 

Please describe anything in the unit that you did not particularly like. 



Please describe anything that you found particularly helpful. (Please use the back of 
this sheet if you need more space.) 



